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ABSTRACT 



Context. The process of high-mass star formation is still shrouded in controversy. Models are still tentative and current observations 
are just beginning to probe the densest inner regions of giant molecular clouds. 

Aims. The study of high-mass star formation requires the observation and analysis of high density gas. This can be achieved by the 
detection of emission from higher rotational transitions of molecules in the sub-millimeter. Here, we studied the high-mass clump 
G30.79 FIR 10 by observing molecular emission in the 345 GHz band. The goal is to understand the gravitational state of this clump, 
considering turbulence and magnetic fields, and to study the kinematics of dense gas. 

Methods. We approached this region by mapping the spatial distribution of HCO + (7 = 4 -» 3), H 13 CO + (7 = 4 -> 3), CS(7 = 7 -» 6), 
12 CO(/ = 3 — > 2), and l3 CO(7 = 3 — > 2) molecular emission by using the ASTE telescope and by observing the 12 C 18 0(y = 3 — > 2), 
HCN(7 = 4 -> 3), and H 13 CN(7 = 4 -» 3) molecular transitions with the APEX telescope. 

Results. Infalling motions were detected and modeled toward this source. A mean infall velocity of 0.5 km s with an infall mass rate 
of 5 X 10~ 3 M Q /yr was obtained. Also, a previously estimated value for the magnetic field strength in the plane of the sky was refined 
to be 855 /iG which we used to calculate a mass-to-magnetic flux ratio, A = 1.9, or super-critical. The virial mass from turbulent 
motions was also calculated finding M vir = 563 Mq, which gives a ratio of M submm /M vir =5.9. Both values strongly suggest that this 
clump must be in a state of gravitational collapse. Additionally, we estimated the HCO + abundance, obtaining X(HCO + )= 2.4xl0~ 10 . 

Key words. ISM: individual objects - ISM: sub-millimeter - Stars: Formation - ISM: Magnetic Fields 



1. Introduction 

It is well known that massive stars (>1OM ) form in giant 
molecular clouds (GMCs), the largest molecular gas and dust 
complexes in our galaxy. The amount of gas involved in this 
process is several orders of magnitude larger than in the low- 
mass star formation scenario, which significantly increases the 
uncertainties that makes its study difficult. In contrast to the low- 
mass star formation case, where theo ry and observations hav e 
established clear evolutionary steps dMcKee & Ostriket1l2007h . 
there are no well defined stages for the evolution of proto- 
massive stars. Presently, two different mechanisms are proposed 
as the dominant process involved in the formation of mas- 
sive sta rs. One mechanism is the coalescence of small frag- 
ments dBonnell et alj Il998t IStahler et all 120001: iBonnell et al] 
The other is accretion directly onto a massive proto- 
stellar object, which is supported b y accumulating evidence for 
disks around proto-mas sive stars (Beuther & Sridharan 2007; 
iMcKee & Ostriketll2007h . 

The rate of formation of massive stars can be up to several 
orders of magnitude lower than low-mass stars, and they ap- 
pear to be born in clusters. The process is highly energetic 
and dynamic, a massive star will quickly perturb and ionize its 
surrounding medium, which makes its study q uite challenging 
(IWood & Churchwelll[T989l: iHoare et al][2007h . Therefore, it is 
crucial to probe the densest regions in high massive star forming 
clumps to understand their physical and chemical conditions. In 



this regard, high rotational levels from molecules emitting in the 
sub-millimeter windows are the ideal tracers to go for. The com- 
bination of optically thick and optically thin isotopomers of the 
same molecular species can be used to infer physical properties 
from star formation sites. 

The magnetic field is likely the least known physical parame- 
ter in star formation. While its presence seems to be ubiquitous 
within the ISM, wit h strengths ranging from a few fiG to mG 
dCrutcher et alJ[l999h . there are surprisingly few observations. 
Moreover, while it is still unclear what role is played by the 
magnetic field in the formation of low-mass stars, the degree 
of uncertainty is even higher for high-mass star formation owing 
to the lack of observations. However, we can say with certainty 
that magnetic fiel ds have been observed towards high-mass sta r 
forming regions dCortes & Crutchetl 120061: iGirart et al.l [2009). 
Therefore, it is important to incorporate it as a relevant physical 
parameter in observations and models of high-mass start forma- 
tion. 

In this paper we present a multi-line study of the high-mass star 
forming region G30.79 FIR 10 The aim of this paper is to un- 
derstand the gravitational state of this clump as well as to study 
the kinematics of the dense gas. We also use previous intefer- 
ometric observations of polarized dust emission to include in- 
formation about the magnetic field toward this region. The pa- 
per is organized as follows, Sect. 1 is the introduction, Sect. 2 
presents the source, Sect. 3 the observational procedure. In Sect. 
4 we present the results, while in Sect. 5 we discuss the abun- 
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dance of HCO + . The kinematics of the gas is studied in Sect. 
5, where we present evidence for infall motions and discuss the 
likelihood of outflows. Additionally, we refined a previous es- 
timation of the magnetic field strength in the plane of the sky 
for this source (Cortes & Crutcherll2006l) . and use it to evaluate 
the gravitational state of this clump. Finally, Sect. 6 presents the 
summary and conclusions. 

2. The source 

G30.79 FIR 10 (hereafter G30.79) is a massive molecular com- 
plex located within the W43 region. It involves an HII region- 
molecular cloud complex near I = 31°, b = 0°, with several 
far-infrared sources, of which G30.79 is the most massive and 
densest component. Figure Q] presents an overview of our ob- 
servations superposed over the 350/mi continuum map from 
iMotte et al] d2003b . lLisz3 ( fT995h observed G30.79 in HCO + and 
13 CO concluding that the presence of several rings and shells 
in the den se molecular gas was a disturbance product of star 
formation. IVallee & Bastienl d2000h observed the dust contin- 
uum emission in this source at 760 yum using JCMT. They found 
linear pol arization of about 1 .9% with a position angle (RA.) 
of 160°. Moo nev et al] d 19951) observed this source at 1.3 mm 
using the IRAM 30 m telescope, detecting a total flux density 
of 13.6 Jy; their wide field map shows the clump and the ex- 
tended HII region in the G30.79 complex. H 2 Q mase rs have 
been observed toward this region (Cesaron i et al. 1988|). which 
are w ithin a half arc-seco nd of the peak in the iMoonev et al] 
d 1995b map. Additionally, lEllingsenl (120071) detected 6.7 GHz 
methanol masers in the proximity of G30.79, which is consid- 
ered to be a signature for massive star formation. No centimeter 
radio-continuum emission seems to be associated with FIR 10, 
suggesting that the source could be in an early stage of evo- 
lution. However, the maser emission already indicat es that star 
forma tion has started and outflows may be present. M otte et al] 
(2003) mapped the W43 main complex in dust continuum emis- 
sion at 1.3 mm and 350/zm with the IRAM 30 m and CSO tele- 
scopes, respectively. They also mapped the HCO + J = 3 — > 2 
line and measured H 13 CO + J — 3 — > 2 towards prominent 
dust maxima. One of the maxima, W43-MM1, corresponds to 
G30.79 and is the compact fragment we observed with ASTE. 
IMotte et al.1 (120031) found Vi sl = QS.Skms -1 , Av = 5.9 kms 1 
(from H 13 CO + ), T dust ~ 19 K, M~ 3600 M , and n(H 2 ) ~ 8.8 x 
10 6 cirT 3 . They estimated the virial mass to be M v i r ~lOOOM , 
suggesting that this compact fragment should be in a state of 
gravitational collapse unless there are other sources o f support 
in addition to kinetic energy. ICortes & Crutcherl (2006) mapped 
G30.79 in dust polarized emission at 1.3 mm using the BIMA 
interferometer founding a polarization pattern, which suggests 
an hour-glass morphology for the field. They also estimated the 
magnetic field strength in the plane of the sky to be ~1.7mG, 
which gave a statistically corrected mass-to-magnetic flux ra- 
tio of 0.9 or critical, where by critical we mean the equilibrium 
value between self-gravitation and magnetic field support for the 
cloud. 



3. Observations 

3.1. ASTE observations 

G30.79 FIR 10 was observed during September 2006 using 
the Atacama Sub-millimeter Telescope Experiment (ASTE) 
from the Nati onal Astronomical Observatory of Japan (NAOJ) 
dKohnol 120 05). The telescope is located at Pampa la bola in 



the Chilean Andes plateau reserve for Astronomical research 
at 4900 meters of altitude. ASTE is a 10 m diameter antenna 
equipped with a 345 GHz double side band SIS-mixer receiver. 
We simultaneously observed 12 CO(7 = 3 -> 2) and HCO + (7 = 
4 3), CS(7 = 7^6) and 13 CO(7 = 3 -> 2), and 
H 13 CO + (7 = 4 — > 3) with a beam size of ~ 22", and a ve- 
locity resolution of 0.1 kms -1 , setting the MAC (which is a XF- 
type digital spectro-correlator) to a bandwidth of 128 MHz. The 
pointing accuracy was in the order of 2", with Uranus used as 
the pointing source. The observations were done by perform- 
ing raster maps (position switching) with a grid spacing of 15" 
(see TableQ] for map sizes), operating the telescope remotely 
from the ASTE base in San Pedro de Atacama under good 
weather conditions (precipitable water vapor or PWV < 1 mm). 
Our reference position was (a, 6) = (18 h 47 m 46.9 s , -1°54'29.1") 
(J 2000), which coincides with the peak du st emission reported 
by Moone v et al] d 19951) : lM~otte et al] (120031) : ICortes & Crutcherl 
(2006). Unfortunately, the off position used to subtract the con- 
tinuum had emission in the I2 CO(7 = 3 — > 2) line. However, 
only the 12 CO line-wings were affected (see discussion below). 
The M17SW molecular complex was used as an intensity cal- 
ibrator. All temperatures are presented as T m b =T^/?7 m b, where 
?7mb = 0.7 1+0.07. Initial data reduction and calibration was done 
using the NEWSTAR package, and the calibrated data were later 
exported into our own software for analysis and plotting. 

3.2. APEX observations 

Observations were performed during the first week of August 
2008 using the Swedish Heterodyne Facility Instrument (SHFI) 
mounted on the Atacama Pathfinder Experiment telescope 
(APEX) dGusten et al]|2006l) . located at llano de Chajnantor in 
the Chilean Andes. We tuned SHFI at 329.3 GHz in order to de- 
tect the I2 C ls O(7 = 3 -> 2) molecular transition, to 354.5 GHz 
for HCN(7 = 4 -> 3), and to 345.3 for H 13 CN(7 = 4 -> 3). 
The spectrometer was set up to 8192 channels with a resolu- 
tion of 0.1 kms -1 . The main beam efficiency is ?/=0.73±0.07 as 
measured by the APEX staff, with a pointing accuracy better 
than 2" and a beam size of 19". Jupiter and R-Aql were used as 
intensity calibrators where the observations were calibrated by 
the usual chopper-wheel method. The observations were done in 
raster mode with spacings of 15" from the same reference posi- 
tion used for the ASTE observations. The initial data reduction 
was done with the GILDAS-CLASS reduction package and the 
final analysis with our own software tools. 

4. Results 

4.1. l2 CO, 13 CO, and 12 C 18 (J = 3^2) results 

The 12 CO (J = 3 — > 2) line was detected all over the sam- 
pled region, covering an area of 3.5' x 2.5' around the peak 
from the dust emission at 1.3 mm (see Fig. 1) The velocity in- 
tegrated emission map is presented in Fig. [2] where the 12 CO 
emission was integrated over the complete available velocity 
range (43 to 154 km s -1 ). The peak in the map corresponds to 
317Kkms _1 located at about 85" to the south-west from the 
reference position, while the peak at the reference position is 
234Kkms~ 1 . Overall we see most of the emission clustered in 
a south-west to north-east elongated pattern, which is consis- 
tent with the morpholog y o f the dust emission se en in the maps 
from lMotte et al.1 (120031) and lMoonev et al.l(ll995l) . However, the 
peak in our 12 CO map does not coincide with the peak in dust 
emission (located at our reference position). Figure [3] shows ve- 
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Fig. 1. Overview of the G30.79FIR 10 massive star forming region. The main panel represents the 350 /mi continuum emission 
from lMotte et al.l (12003). Figure 1. Contours and grayscale are 12, 24, 48, 96, and 200 JybeairT 1 . The cluster of massive stars is 
indicated by the star symbol, while masers are indicated by the symbols presented at the upper right of the figure. Superimposed are 



our new 12 CO(7 = 3 — > 2) and CS(J = 7 — > 6) integrated velocity maps to show the extent of our observations. 



Table 1. Line parameters 



Line 


Transition 


Frequency 


Observed Area 




J Tdv 


Vpeak 


AV 








[GHz] 


[arcsec 2 ] 


[K] 


[Kkms- 1 ] 


[kms- 1 ] 


[kms" 1 ] 


HCO + 


(/ = 


4^3) 


356.7342880 


200" x 150" 


5.7 ±0.1 


42.6 ±1.8 


94.8 


8.4 


H 13 CO + 


u = 


4^3) 


346.9983381 


90" x 90" 


1.2 +0.02 


5.5 ±0.6 


98.1 


3.0 


HCN 


(/ = 


4->3) 


354.5054779 


200" x 150" 


3.8+0.1 


73.5 ±13.8 


93.8 


9.5 


H 13 CN 


(/= 


4^3) 


345.3397750 


30" x 30" 


1.4 ±0.1 


14.8 ±0.4 


99.1 


4.8 


cs 


u = 


7^6) 


342.8828503 


90" x 90" 


3.6 ±0.2 


38.7 ±0.9 


98.5 


3.8 


12 C 16Q 


(J = 


3^2) 


345.7959899 


200" x 150" 


16 ±2.2 


317 ±21.1 


98.8 


9.6 


13 C 16 


(/= 


3^2) 


330.5879652 


90" x 90" 


8.5 ±1.2 


125 ±7.6 


97.4 


6.3 


12 C 18 


(J = 


3->2) 


329.3305525 


30" x 30" 


3.1 ±0.2 


42.9 ±10.0 


98.8 


4.7 



Parameters for all molecular line observations. The values for T pcak and AV are calculated from Gaussian fits to the central pointing or 
(18 h 47 ra 46.9 s , -1°54'29.1"). For HCN and HCO + , the V pea k was taken from Gaussian fits done to the most intense component, while the 
line-width corresponds to the whole line. For 12 CO and 13 CO the fit was done over the whole spectrum. 



locity channel maps binned every 5 km s _1 (or 50 channels) with 
a resolution of 5.4km s _1 for each map with respect to the orig- 
inal resolution of O.lkms 1 . The channel maps are presented 
in increasing velocity from 84 to 106kms~ 1 covering the most 
intense features in the 12 CO emission. By examining Fig. [3] a 
velocity gradient from the south-west to north-east can be seen 
with each successive channel where the peaks in emission are 
at 89 and 95kms _1 and at almost opposite spatial locations. 
Particularly interesting is the inten se emission seen at the south- 
west part of the map. According to lMotte et al.l (I2003L see Fig. 3 
and also our Fig. Q]), the giant HII region produced by a Wolf- 



Rayet cluster of ma ssive stars has not ye t reached the G30.79 
clump. However, the lMoonev et alj d 1995b maps put the HII re- 
gion at an interface with the clump. Because we could not sam- 
ple beyond 105" west, we can only speculate about the nature of 
this intense 12 CO(7 = 3 — > 2) emission and whether it is related 
to an interface between this clump and the H II region without 
further observations. 

To study the velocity structure at selected locations in the in- 
tegrated emission map, we present 12 CO spectra from three posi- 
tions in the map in Fig. [4] These positions correspond to averages 
over nine pointings centered at (Aa, A6) = (0, 0), (-90, -40) and 
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Fig. 2. Velocity integrated 12 CO(7 = 3 — > 2) map. The axis units are offsets in arc-seconds from the reference position 
(18 h 47 m 46.9 s ,-l°54'29.1") in J2000 coordinates. The gray-scale is in Kkms -1 , while the contours are 3, 4, 5, 6, 7, 8, 10, 12, 
and 14<x with cr — 20 Kkms -1 . 
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Fig. 3. Channel maps for the 12 CO line emission. Maps are shown as averages of 50 channels giving a velocity resolution of 
5.4km s~'for each map. The average velocity in each plot is shown by the box in lower left corner of each map. Both Aa and A6 
are arc-second offset with respect to the reference position. The pixel scale is in K degrees and is normalized to the maximum and 
minimum values from all the maps set. The contours are from 3, 6, 9, 12, and 15cr, where cr = 0.9 K. for all maps and calculated 
from a velocity range with no emission. 



(30,40) arcsec. The (0,0) offset plot corresponds to the spec- 
trum over the reference position and had superposed 13 CO and 
12 C ls O spectra as well. The emission is strong with a peak of 
16 K for the 12 CO, 8.5 K for 13 CO, and 3.1 K for 12 C lg O, ob- 
tained by fitting Gaussian profiles. The remaining spectra cor- 
respond only to 12 CO and are taken from the maxima seen at 
89 and 95 km s _1 velocity channel maps. We did not obtain data 
from the other CO isotopomers at these positions. We immedi- 
ately appreciate that the 12 CO line profiles are complex and have 
many components with the emission appearing to be optically 
thick. But owing to the large angular extension of the whole 
W43 region, we did not find a suitable off position in our ob- 
servations. This is clearly seen in the dip features at the line 
wings of our spectra (see Fig. [4] upper panel). While the emis- 
sion at the off position is weak compared to the peak of the line, 
it is strong enough to disrupt the line-wings. In principle this 
does not affect the integrated intensity map, but it will prevent 
us from concluding about possible outflow motions from the 
CO line-wings. We also see that all spectra show blue-shifted 
peaks relative to Vi sr =98.8 kms -1 . The associated peak veloci- 
ties are 94.24km s" 1 for (0,0), 89.6910ns- 1 for (-90,-40), and 
95.00 km s -1 for (30,40). All spectra show broad line-wings, par- 
ticularly at offset (-90,-40), and a high velocity component at 
v= 1 1 5 km s~ 1 . This peculiar emission may indicate outflow mo- 



tions, but as mentioned before, it is difficult to be decisive due to 
contamination from the off position. 

In mapping 13 CO {J — 3 — > 2), we were only able to cover 
an area of 1.5' x 1.5' around the reference position. The max- 
imum is found at 124.7 Kkms -1 located at 30" west from the 
center. The emission also seems to follow a south-west to north- 
east orientation like 12 CO do. Yet, we cannot confirm the cor- 
relation due the lack of equal coverage for our 13 CO. Fig. [4] 
first panel, shows a spectrum corresponding to the nine central 
pointings, or 30" x 30", in 12 C ls O and 13 CO overlaid on the 
12 CO spectrum. The same complex features seen in 12 CO are 
seen in the 13 CO line profile, particularly the high velocity com- 
ponent at 1 15 km s -1 . Signs of self-absorption near the center of 
the line seem to be present, but not the strong absorption seen in 
the 12 CO profile at the edges of the line wings. From these fea- 
tures, it is likely that the 13 CO emission that we detected towards 
the center is also optically thick. In contrast, 12 C ls O is single- 
peaked, with an almost Gaussian profile, showing no trace of 
self-absorption, which suggests a likely optically thin emission. 

An additional feature of the CO data is a consistent blue- 
shifted peak emission. It has been suggested that optically thick 
lines with a blue-shift ed profile may indicate infall motions 
(lLeung&Brownlll977l) . Both 12 CO and 13 CO emission are sim- 
ilar in complexity showing the same gradient in velocity from 
south-west to north-east, broad line wings, and blue asymme- 
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Fig. 4. Spectra from 12 CO (J = 3 — > 2) are presented. Each 
spectrum corresponds to an average of nine pointings centered 
at offsets, from upper to bottom, (Aa, A£)=(0,0), (-90,-40), and 
(30,40), where the offsets are in arc-seconds. The velocity reso- 
lution for all the spectra is taken to be 0. 1 km s _1 . The first panel 
also shows the spectra from 13 CO (J = 3 — > 2) in red and C ls O 
(J = 3 — > 2) in green. Clearly visible is the negative dip feature, 
likely due to emission in the off position, seen in both the upper 
and lower panels for the 12 CO (J — 3 — > 2) spectra. The compo- 
nent at 1 15 km s _1 is also seen in both 12 CO and 13 CO but not in 

c 18 o. 

try in the line profile (Vi sr =98.8 km s~'). Particularly interest- 
ing is that the line profiles in addition to their blue-shifted peak 
are both s elf-absorbed, which is a lso often seen in infall can- 
didates dKlaassen & Wilso"nll2007h . Considering that our offset 
(0,0) corresponds to the peak in dust emission at both 1.3 mm 
and 350yum, which covers most of the continuum spectrum 
from pre-stellar cores, it is likely that the condensation we are 
seeing is actively forming stars. Unfortunately, we cannot re- 
solve the number of fragments due to the coarse 22" beam of 
the ASTE telescope. Our APEX observations, with their 19" 
resolution will not help either because previous interferomet- 
ric, 4" resolution, observ ations also did not resolve the core 
(ICortes & Crutcherl l2006h . However and due to the expected 
multiplicity in high-mass star forming cores, it is unlikely that 
we are seeing only a single object. 

4.2. HCO + (J = 4^3) and H l3 CO + (J = 4 -> 3) results 

We observed HCO + (7 = 4 -> 3) simultaneously with 12 CO 
(J = 3 — > 2), which also covered an area of 3.5' x 2.5'. The 
velocity integrated emission map is presented in Fig. [5] The mor- 
phology of the emission follows the same orientation south-west 
to north-east as the 12 CO observations, but without the extension 
of 12 CO. Clearly, the strongest emission is clustered at the center 
of the map, to the north-east, while some hint of emission over 
3<r seems to appear at the south-west, roughly at {Aa, A£)=(- 
80", -40"), where the strongest 12 CO emission is located. The 
map m orphology is also co nsistent with the HCO + (/ = 3 — > 2) 
map of iMotte et alJ d2003l see Fig. 2). Figure [6j central panel, 
shows the spectrum from the reference position for both HCO + 
and H 13 CO + . The line profile is clearly non-Gaussian, showing 
evident self-absorption. It is likely that the HCO + emission is 
optically thick. No hint of emission is seen at the off position, 
which is same off position used for 12 CO. Superposed in thick 




80 100 120 140 



Velocity [km/s] 

Fig. 6. High density molecular gas spectra from the G30.79 FIR 
10 clump. The upper panel presents the HCN(7 = 4 — > 3) spec- 
trum. In thick red lines the H 13 CN(7 = 4 — > 3) spectrum is 
superposed. The central panel presents the HCO + (7 = 4 — > 3) 
spectrum. In thick red lines the H 13 CO + (7 = 4 — > 3) spectrum 
is superposed. The peak brightness temperature of H 13 CO + is 
1.2 K. The lower panel presents the CS(7 = 7 — > 6) spectrum. 
All spectra were taken at the reference position. Both HCO + and 
the HCN spectra were binned every five channels, giving a ve- 
locity resolution of 0.54 km s -1 , while CS, H 13 CO + and H 13 CN 
were binned every ten channels, giving a velocity resolution of 
1.1 kms -1 . 

lines is the H 13 CO + (/ = 4 — > 3) spectrum with a peak brightness 
temperature of 1.2K at V=98.1kms -1 . The H 13 CO + emission 
is only significant (over the 3<x level) at the reference position. 
The velocity associated with the peak can be considered to be 
at the Visr within the boundaries of the bin. The emission is also 
coincident with the strong absorption dip in HCO + , suggesting 
self-absorption. Additionally, because the line profile appears to 
be Gaussian in shape, it is likely that the H I3 CO + emission is 
optically thin. 

4.3. CS(J = 7^6) results 

The integrated velocity map for CS(7 = 7 — > 6) is presented in 
Fig- El As for H 13 CO + , the CS emission is fairly compact, aris- 
ing only from the center of the map, which suggests that most 
of the activity is coincident with the peak of the dust emission. 
The spectra from the central 30" x 30" are shown in Fig.|6]in the 
lower panel. The CS emission presents a peak brightness tem- 
perature of 3.6K at V=98.1 kms^'with some excess emission 
in its line-wings, which may be due to outflowing motions. Note 
that the CS emission is not self-absorbed as seen with HCO + and 
HCN. It is possible that the molecule is not abundant enough to 
become self-absorbed. However, another possibility is that the 
CS gas is bound inside the dense core, which would explain why 
its emission is not widespread (as indicated by our map). 

4.4. HCN(J = 4^3) and H l3 CN(J = 4^3) results 

The HCN(7 = 4 — > 3) was observed in position-switching 
mode over and area of 200" x 150" to sample the same re- 
gion as HCO + . We looked for emission only in the most sig- 
nificant areas within the region. The H I3 CN(7 = 4 — > 3) was 
only observed over the central 30" x 30", where we expected 
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Fig. 5. Velocity integrated emission map for HCO + (7 = 4 — > 3) over the whole velocity range (45 to 150 km s ') is shown. The 
pixel scale is in units of Kkms" 1 , while the contours are on levels of 3, 4, 6, 8, and 10 x cr, with cr = 3.7 Kkms -1 . 
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Fig. 7. Velocity integrated emission map for CS(7 = 7 — > 6) over the whole velocity range (45 to 150km s '). The pixel map has 
units of K kms -1 . The contour levels are 3, 6, 9, and 15 xc, with cr = 1.7 Kkms" 1 . 



to find most of the emission. Figure [8] presents the velocity in- 
tegrated emission map for HCN. While the emission is mostly 
concentrated at the center, there is a hint for a gradient along 
the south-west to north east direction as with HCO + and CO. 
The extent of the emission appears to be midway between the 
HCO + (7 = 4^3) and the CS(7 = 7^6) where the CS map 
presents the most compact morphology. Fig.|6]shows in its upper 
panel the HCN(7 = 4 -> 3) overlaid by the H I3 CN(/ = 4 -> 3) 
spectrum in red. The HCN emission is clearly opticaly thick 
with a self-absorption feature as indicated by the optically thin 
H 13 CN. The double peaked HCN spectrum has a stronger blue 
peak at 93.8 kms" 1 ; while the H 13 CN is peaked at 98.3 kms" 1 , 
both taken from Gaussian fits. This type of spectrum is often 
seen toward infall candidates like the well-studied source B335 
dZhou et al 1 119931) . 

5. Analysis and discussion 

5.1. The HCO + abundance 

In the current paradigm, it is suggested that carbon bearing 
molecules will freeze onto grai n surfaces, forming ices, in 
cold and dense cores (see e.g. Bergin & Tafalla 2007); while 
some other molecular abundances, like HCO + , might get en- 
hanced by non-therm al motions such as outflows and/or infall 
dRawlings et alJl2004l) . In order to investigate the abundance of 



HCO + toward this region, we estimate the physical parameters 
associated with this molecule. The optical depth for HCO + is 
estimated by the ratio between the opt ically thick and optically 
thin brightness temperatures by using (IChoi et al.ll 19931) 



l mb4^3 
r mb(13)4->3 



1 



I _ e -U 3 /X ' 



(1) 



where r m b4^3 is the brightness temperature from the HCO + line, 
?mb(i3)4->3 is the H 13 CO + line temperature, and X is the ratio be- 
tween X=[ 12 C]/[ 13 C] taken to be 50, a value r epresentative for 
the m olecular ring at the 4 kpc inner galaxy (Wil son & Rood! 
1994). Because the HCO + profile is self-absorbed, there is a sig- 
nificant uncertainty in deriving the line intensity. Two possible 
approaches can be followed: either we pick the strongest peak, 
or a Gaussian profile can be adjus ted to the whole lin e by mask- 
ing out the self absorption dip. iPurcell et al.1 (12006) found no 
major differences from these two approaches in their results for 
a sample of 40 high-mass star forming regions. Thus, we used 
the strongest peak for the HCO + line intensity in this calcula- 
tion. From these numbers we obtained tJ| 3 o+ = 11.3 indicating 
an optically thick HCO + line, which is to be expected due to the 
self-absorption seen in the line profile (see Fig. |6j. Also, note 
that both the HCO + and the H 13 CO + line were obtained from 
the reference position and with a beam size of about 18", which 
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Fig. 8. Velocity integrated emission map for HCN(7 — 4 — > 3). The integration is done over 60 to 140 km s~'to include all the broad 
HCN emission. The pixel map has units of K kms -1 . The contour levels are 3, 6, 9, and 15 xcr, with <x = 3.0 K kms -1 . The color 
scale indicates the emission level with white used for positions not observed or with zero emission. 



encloses the dust cor e detected from the interf erometric obser- 
vations of this region dCortes & Crutcherif 2006. see Figure 1). 
To o btain the column density fo r the upper rotational level, we 
use (iGoldsmith & Langerll 1999b 

where v is the line frequency, A43 is the Einstein spontaneous 
emission coefficient (3.6269 x 10 -3 s -1 for HCO + (7 = 4 -> 3)), 
Af2 fl is the antenna solid angle and A£2 S is the source solid angle 
for the core which gives a beam filling factor of 0.73 for HCO + . 
The total column density is then given by 



where Z is the partition function, g — 2J + 1 is the statistical 
weig ht of the upper level, E u = 42.8 K is the energy for the upper 
level (ISchoier et al.ll2005l). and T = 30 K is the dust temperature 
used bv lCortes & Crutcher (2006) to calculate the total hydrogen 
column density NQI2). This is justified under the LTE assump- 
tion. For linear molecules, the partition function is well approx- 
imated by Z = kT/hBo, whe re Bp _ 45 GHz is t he rotational 
constant for HCO + taken from lDefrees et al](11982l) . In this way, 
we estimated a HCO + column density N(HCO + ) = 1.5 xlO 14 
cm -2 . From a sample of 40 high-mass star forming region s se- 
lected from methanol maser detections. iPurcell et al.l f2006) esti- 
mated the HCO + column densities ranging from 2.5xl0 14 cm -2 
to 71 XlO 14 cm~ 2 . Our column density estimation is consistent 
with the smallest values obtained in their work. Now, by using 
the column density of hydrogen derived from dust emission or 
NQI2) = 6.5 xlO 23 cirT 2 , we calculate a total abundance for 
HCO + or X(HCO + ) = 2.4 xlO -10 . As previously mentioned, it is 
expected that carbon bearing molecules will get depleted toward 
cold and dense regions. We can quantify this by calculating the 
depletion factor fo defined as the ratio betw een the average frac- 
tional abundance and the observed one. iLucas & Lisz tl d 19961) 
and [Liszt & Lucasl (120001) determined a Xism =2x10 9 from a 
sampling of HCO + (7 = 1 -> 0) toward the diffuse ISM in ab- 
sorption against background QSOs. By using this value for Xism, 
we find fo = 8.4, which suggests that a moderate depletion of 
HCO + toward this clump when compared to othe r high-mass star 
forming regions such as G305 with fo = 1.5 ( Wals h & Burtonl 
2006; Ber gin et alJll9 97). However, we maybe underestimating 
the HCO + column density due to the uncertainty introduced by 
the brightness temperature used from the observed self absorbed 



profile, which would lead a lower depletion value for HCO + . 
These results, along with the lack of infrared sources detected at 
this clump, and with no compelling evidence for outflow emis- 
sion, even though maser emission has been observed (see next 
section), suggest an early (pre-hot core) evolutionary stage for 
G30.79 FIR 10. In an advanced stage, it is expected that the ra- 
diation coming from the central condensation will evaporate the 
ice mantles increasing the abundance of molecules like HCO + , 
creating what is known as a hot core. Note though, that the de- 
tection of maser emission toward this clump suggests that the 
star formation process has already started. 

5.2. An outflow in G30.79 FIR 10? 

One of the signatures of high-mass star formation are the 
powerful molecular outflows observed t oward these regions 
(IShepherd et~aD I2007I: iBourke et alJll997l) . Because FIR 10 is 
a massive core, it is likely that outflowing motions are present 
or will develop over time. However, to untangle the outflow mo- 
tions from our molecular emission observations is certainly not 
trivial. The choice of sub-millimeter emission lines as a tool 
to study this region allow us to separate the most dense com- 
ponents from the rest of the molecular core. Particularly, 12 CO 
(J — 3 — > 2) has been successf ully used to trace outflowing mo- 
tion from star forming regions (IChoi et al.lll993l However, our 
I2 CO (J — 3 — > 2) results are inconclusive because of emission 
at the off position, which affected the line wings in our spec- 
tra. Also, the high velocity component seen at V= 115 kms -1 , 
which is ~20kms _1 over the V] sr and is only seen in 12 CO and 
13 CO but not observed either in HCO + , HCN, or CS. It is likely 
that this component is not dense enough to excite either of these 
lines; indeed, this high velocity component might be optically 
thin. Although it could be interpreted as high velocity outflow 
emission, its widespread spatial distribution makes this interpre- 
tation un-likely. Another possibility is that this component corre- 
sponds to a foreground object, or cloud, which is not associated 
with this clump. Additionally, it is difficult to accelerate the gas 
to such high velocities without dissociation. 

Outflows are discovered through their signature in the line- 
wings of molecular emission lines. Fig.|9]shows the HCO + and 
CS spectra with their corresponding Gaussian fittings and the 
residual emission in their line-wings. Between 100 and 110 
kms we found some significant emission over the 3<x level, 
with cr = 0.08 K, in both CS and HCO + spectrum. However, 
the blue-shifted part of the emission, for velocities lower than 
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95 km s , does not show significant traces of residual emission. 
Even though the red-shifted excess emission may be due to non- 
thermal motions, the bipolar nature of an outflow must be clearly 
stated, which we cannot do with these data. Therefore, we can- 
not yet come to a conclusion about the presence of a bipolar 
outflow towards this source. Our lack of spatial resolution due to 
the distance to G30.79 FIR 10 might be the reason behind this. 
However, outflows are ubiquitous in the high-mass star form- 
ing regions, so it is likely that they are present or will develop 
over time. Additional observations of shocked excited chemistry, 
such as SiO, SO, or SO2 may help confirming an outflow in this 
clump. 

5.3. Infall motions 

It is not clear whether high-mass stars forms through accretion 
or through a different process such as coalescence of less mas- 
sive fragments. This situation is difficult to distinguish due to 
the physical complexities involved in the evolution of high den- 
sity gas and dust. The study of the kinematics and dynamics at 
the earliest phases, along with the detection of accretion disks, 
could clarify this uncertainty. In this scope, the determination of 
infalling motions is a first step toward the identification of col- 
lapsing pre-stellar objects. The characterization of these motions 
is a challenging topic, the current avenue towards investigating 
infalling is the study of asymmetries in molecular line profiles. 
iLeung & Brownl dl977h suggested that an asymmetry in the line 
profile toward the blue may indicate the presence of infalling 
motions. Thus, the low excitation, red-shifted infalling layers of 
gas in the front part of the cloud absorbs some of the emission 
from the rest of the gas. This red-shifted self-absorption is what 
makes the spectrum present a brighter blue peak. To quantify 
this asymmetry, the calculation of the norm alized line velocity 
difference has been used by many a uthors dFuller et al.l 12005; 
ISzvmczak et al.ll2007tlWuetaIll2007h . 



Table 2. Parameters from the model fit 



AVbe 



T/max _ T/max 
thick thin 

AVthin 



(4) 



where V™* k is the velocity from the peak emission of the op- 
tically thick line, V t '™ n x is the velocity from the peak emission 
of the optically thin line, and AV t hi n is the FWHM line width 
from the optically thin molecular tracer. A negative AV^e will 
correspond to blue-shifted gas velocities, which are indicative of 
infalling motions. We calculated AVt> e for the high density trac- 
ers HCN, CS, and HCO + using H 13 CO + and H 13 CN as the op- 
tically thin tracers. These molecules have high dipole moments 
and therefore require higher critical densities to get excited. In 
particular we used H 13 CO + as the thin component for HCO + 
and H 13 CN for CS and HCN to keep ion-molecules and neutral 
molecules separate. Table[3]gives the values of AVb e for all three 
lines. The difficulty in separating infalling evidence from the line 
asymmetry arises because other dynamic phenomena such as ro- 
tation and outflows will also produce red and blue asymmetries 
in the molecular profiles. However, infalling is the only motion 
which will produce only blue asymmetries. The expected infall 
signature will manifest in spectral lines with a double peaked 
profile where the blue peak is stronger than the red peak and 
the dip is due to self absorption. A blue asymmetry is seen in all 
spectra, with the highest value in HCO + , but the characteristic in- 
falling double peak profile is only seen in both HCN and HCO + . 
Both appear to be self absorbed and have a stronger blue compo- 
nent (see Fig. [6]). The case of HCO + is interesting, the red com- 
ponent is almost non-existent having the largest blue asymmetry 
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The parameters f rom the best fit obtained by adjusting the 
iDi Francesco et alj d2001l) infalling model to our HCN(7 = 4^3) 
spectrum. 



in all three molecules. The CS spectrum is single peaked to the 
blue, but does not show the self-absorbed features. It has the low- 
est intensity, which could be cause by carbon and sulfur deple- 
tion in this obje ct. In this way , it may not have enough abundance 
to self-absorb. lEvansI d2003l) suggested a procedure to evaluate 
whether infall is likely or not in a star forming core. We followed 
the procedure checklist and found that our HCN observations, 
taken at the reference position, meet the conditions required for 
infall. Moreover, the same properties are observed in the HCN 
spectra from the surrounding pointings (see Fig.fTOii. Almost all 
of them present double peaked profiles with a stronger blue peak, 
clearly suggesting infalling. Together with the previous checklist 
we can conclude that our HCN spectra meet the requirements to 
suggest infall movements toward this source. Thus, we quanti- 
fied the infall parameters through modeling . To do this, we ap- 



plied the simple infall ing model devised b y M yers et"afl (fl996) 
and later improved by Di Frances co et al.l (12001). In this model, 



the clump is approximated by two infalling gas layers, a front 
and a rear layer, with a central source modeling the pre-stellar 
core. Thus, the observed brightness temperature is quantified by 
the following equation, where the subscripts "f ', "r", "C", and 
"cmb" stand for the front layer, the rear layer, the central source, 
and the cosmic background. 



AJ B = (Jf - /c)[l - e- T '] + (1 - 0)(7 r - 7 cmb )[l - e 



<T,+T,) 



]■ (5) 



The main terms in this equation are the Planck excitation tem- 
perature given by 7; = To /[exp {To /T{) - 1] with To = hv/k, 
and T\ corresponding to either Tf, T T , T c , and r cm b- Also, 7c = 
07 c + (1 - <t>)y r where <L> is the beam filling factor of the con- 
tinuum source, which due to our large beam size is assumed to 
be 0. The t; expressions correspond to the optical depths , which 
we assumed to be Gaussian (following Mye rs et alJ 19961) . Thus, 
the front and rear optical depths are given by 



T f = t exp 
T r = T exp 



-(V - Vf - v lsr ) 2 



2cr 2 

-(V + V r - Vi sr ) 2 



2cr 2 



(6) 
(7) 



where to is the peak optical depth for both the front and the rear 
layers, Vf and v r are the infalling velocity for both slabs and cr is 
the velocity dispersion. The model was fitted in a two step pro- 
cess, first by doing a multi Gaussian fit component, which pro- 
vided the input parameters to later adjust Eqn.[5]to our data by 
minimizin g the y 2 fun ction through the Levenburg-Marquardt 
algorithm (Press 2002). The best fit to the HCN spectrum is 
shown in Fig. [TTJ and the fit parameters are presented in Table 
5.31 Note that this is the simplest possible model that we can use 
to characterize infalling gas. Moreover, the structure of the HCN 
line is complex and, even though the features seen in the line 
wings may come from outflowing motions, we did not consider 
them in this simple model. 
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Fig. 9. Residual emission after Gaussian substraction is shown in both panels. The lower panel shows the HCO + (7 = 4 — > 3) 
profile with an adjusted Gaussian to the left. To the right, in the same panel, the residual emission from the Gaussian subtraction is 
presented. The upper panel shows CS(7 = 7 — > 6) profile with the adjusted Gaussian profile and its residual emission to the right. 
The cr value is calculated to be 0.08 K. Both spectra were taken as averages over an area of 30" x 30". 
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Fig. 10. The central 48" x 48" region compose of 9 pointings separated by 15" with a beam size of 18" each is shown. The green 
line represents the systemic velocity at 98.8 km s _I ; while the scale is given in K (main beam temperature) 




Fig. 11. Central HCN spectrum from G30.79 FIR 10. Overlaid in thick lines is the best fit for the infalling model. 
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Table 3. Normalized velocity difference. 



Line 


AV be 


HCN(7 = 4 -> 3) 


-1.1 


HCO + (y = 4^3) 


-1.2 


CS(J = 7^6) 


-0.1 



The asymmetry of line profiles is presented by the calculation of 
AV bc , as shown in Eqn. [4] where the optically thin species correspond 
to H 13 CO + for HCO + and H 13 CN for HCN and CS. The values 
used to compute the blue excess are taken from Gaussian fits to the 
corresponding spectra. 



The infalling velocity for both slabs is estimated by the infalling 
model to be Vf = 0.7 km s 1 and v r = -0.3 kms" 1 . A simple 
estimation can also be done by using the analytical model de- 
rived by iMvers et"afl ([1996) for a contracting cloud. By using 
Eqn. 9, from their model, we estimated an infalling velocity of 
0.5 kms 1 . This result is in good agreement with the the val- 
ues obtained by adjusting the model to our data. The next step 
is to estimate the mass i nfall rate, which we do by following 
dKlaassen & Wilsonll2007h 

AM M pVv'i n 4 , 
M = — — ~ — = — — = -7rn H2 fm H R l v in , (8) 
at t R 3 

where fi = 2.35 is the mean molecular weight, R = 0.3 pc 
is the geometric radius taken from the interferometric obser- 
vations, and «h, is the gas density. The infall rate obtained for 
this source is about 5 x 10~ 3 M /yr. Both infall velocity and 
infall mass rate are consistent with the values seen at other 
high-mass star forming regions, like G 10.47, which has an es- 
timated irrfidl_velocit^_ofM of 10~ 2 
M Q /yr dKlaassen & Wilsonl2007l) . lBarnes eF al. (2008) estimated 
infall rates of 3 x 10~ 2 M /yr with an infalling velocity of 
0.3 kms^'forthe G286.21+0. 17 molecular clump. Both regions 
present similar gas densities as G30.79 FIR 10 (1.9 x 10 5 cm -3 
for G286.21+0.17 and 7.2 x 10 s cm" 3 for G 10.47) and similar 
mass estimations (4.1 x 10 4 Mpf or G286.21+0 17 an d 7.3 x 10 3 
M Q for G10.47 as estimated by IShirlev et"ai]|2003l) . These re- 
sults are consistent with our findings where these regions appear 
to be in a similarly evolutionary stage. 

5.4. Dynamical state of G30. 79 
5.4.1. The magnetic field 

The dynamical state of high-mass star forming cores is difficult 
to establish. The accumulation of a large amount of molecular 
gas is still poorly understood, as is what triggers the gravita- 
tional collapse. Moreover, it is known that magnetic fields are 
present in these regions with stre ngths substantially larger than 
in low -mass star forming regions dLaill200lt ICortes & Crutcherl 
2006). Thus, the magnetic field can have a profound effect on 
the dynamical evol ution of the gas as postulate d by ambipolar 
diffusion theories dCiolek & Mouschovias| [T993). Our object of 
study is not an exception to the above. Magnetic fields have been 
observed through polarized emission from dust to wards this core 
dVallee & BastienfeOOrilCortes & Crutcherl2006l) . However, it is 
not yet clear how important the field is with respect to turbulence 
in the high-mass star forming process, but its presence alone is a 
reason to consider its effects. 

From the interferometric obse rvations of polarized dust 
emission, ICortes & Crutcherl d2006l) found a smooth polarimet- 
ric pattern centered at the G30.79 FIR 10 clump. Assuming mag- 



netic alignment of dust grains, they proposed an hourglass mor- 
phology for the field. Additionally, they estimated the magnetic 
field strength o n the plane of the sky using th e Chandrasekhar- 
Fermi method (Chandrasekhar & FerrnJ[l953|), finding a value 
of B pos = 1.7 mG. In this way, the mass-to-magnetic flux ratio 
was calculated to be A = 0.9 or critical. In this work, we refined 
the previously described estimation by considering our molecu- 
lar line observations. Before stating this improvement, we will 
briefly introduce the Chandrasekhar-Fermi method by the fol- 
lowing equation. 

Spos = 9.3 \ 2 ' , (9) 

Of 

where B pos is the magnetic field strength in fiG, n(H2) is the vol- 
umetric gas number density in cm 4 , AV is the FWHM velocity 
width in km s" 1 , and 54> is the polarization position angle disper- 
sion in degrees. lCortes & Crutcherl d2006l) used n(H 2 ) = 4.5xl0 5 
cirT 3 calculated from their dust emission, while 5(f> = 21°. 9 was 
obtained from the polarization data. The line-width, AV, was 
taken f rom the H 13 CO + (7 = 3 -> 2) observations bv lMotte et al.1 
d2003l) for the G30.79 clump (called MM1 and listed at Table 2 
in their work). It is here where we improve the estimation by 
using AV derived from our H 13 CO + (7 = 4 — > 3) observations. 
The (J = 4 — > 3) transition traces higher density molecular gas, 
while being better coupled to the field than H I3 CN or other neu- 
tral species. Also, the H 13 CO + (7 = 4 — > 3) emission is optically 
thin, which might not be the case for the {J — 3 — > 2) transi- 
tion. Additionally, this higher transition comes from a smaller 
beam size, which engulfs only the FIR 10 clump (of about 20" 
in size), while the (J — 3 — > 2) transition comes a higher 
beam size (~ 30"), which might be including velocity compo- 
nents not associated to the field perturbations. Thus, by using 
AV = 3.0kms _1 we obtain a value for Z? pos =855 /iG, which 
gives a statistically corrected mass-to-magnetic flux of A = 1 .9, 
or super-critical, consistent with the infalling scenario suggested 
in this work. Thus, the magnetic field is not strong enough to 
support the core. 

Whil e the polar ization pattern observed by 
ICortes & Crutcherl d2006l) is smooth, it is unclear how dis- 
torted by turbulence the field might be. This source is 5.5 
kpc away which gives 0.03 pc/1" or 0.1 pc/beamsize. This 
length scale is larger than the predicted turbul ence length 
scale of 1 mpc dYan et alJl2004l iLi & Houdell2008l) . Therefore, 
the dispersion in the field lines, or the turbulent component 
of the field, cannot be directl y seen from the polarization 
pattern. Hildebrand et al.l J2009) developed an approximation 
to calculate the plane of the sky dispersion function for the 
magnetic field in turbulent molecular clouds. By using polarized 
emission from dust, they calculated a dispersion function for the 
polarization position angle as a quadratic function of the length 
scale. Two main assumptions are made in this calculation. 
The length scale is larger than the turbulence, or correlation 
length scale, and smaller than the length scale associated with 
large scale variations of the field. Our BIMA polarization data 
satisfy both assumptions, the interferometric observations have 
a length scale on the order of 0.1 pc while the field varies 
over 1 pc distances (as seen from the polari zation pattern). 
This approximation is stated in Eqn. 2 in Hildebrand et al.l 
(2009). From this approximation, the ratio between the turbulent 
component of the field on the plane of the sky and the main 
magnetic field strength can be estimated as 
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Fig. 12. The polarization position angle dispersion function is shown here in closed circles along with its respective error bars. The 
dispersion function points are join with thick lines, while superposed is the fit for EquationfTOlto the first 3 points. 
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where b is the turbulent contribution to the dispersion function, 
{b^i is the turbulent component of the field, and Bo is the main 
magnetic field. We applied this method to our BIMA data by cal- 
culating the dispersion function for polarization position angle as 
shown in Fig. [12] The value of b — 0.42 is given by the intersect 
of the best fit to the data points, and gives a turbulent-to-main 
magnetic field ratio of {b^ /Bq = 0.3. This result suggest that 
30% of the total magnetic field strength is in the turbulence. 



we obtained the ratio of the virial mass to the dust mass to be 
M su bmm/Mvir = 5.9. Values of Msubmm/Mvir > 0.5 are consid- 
ered t o remain gravitationally bound according to lPound& Blitz 
(119931) . At the same time, we recall our previous calculation of 
the mass-to-magnetic flux ratio, which was found to be super- 
critical. This means that there is not enough energy in the turbu- 
lence and/or in the magnetic field to balance gravity. Therefore 
and even though the uncertainties involved in the total mass esti- 
mations for this object, the results suggest that G30.79 FIR 10 is 
bound and must be undergoing gravitational collapse. The whole 
picture appears to be self-consistent, we have detected infalling 
motions into the clump where neither turbulence or the magnetic 
field have enough energy to counterbalance gravity. 



5.4.2. The gravitational equilibrium in G30.79 FIR 10 

This source is one of the few where information about the mag- 
netic field and turbulence is present. Thus, we should be able 
to characterize with confidence the state of gravitational equi- 
librium in this source. In orde r to quantify this, we wil l fol- 
low the analysis as state d by iMcKee & Zweibeil dl992l) and 
iBertoldi & McKed (I1992I) : 

2(7" - To) + A4 + W = (11) 

where T is the total kinetic energy including thermal and non- 
thermal motions, To correspond to the surface term for the ki- 
netic energy, or external pressure, M is the magnetic energy, and 
W is t he gravitational energy term. We will follow iMotte et aU 
(2003) and assume the giant HII region present in the W43 com- 
plex does not reach the clump. In this way, we can drop To 
out the equation. Only the kinetic and magnetic terms remains 
to balance gravity and as was said before, we have information 
about both of them. 

Motte ^tall (120031) estimated the virial mass without know- 
ing the magnetic field, stating that G30.79 FIR 10 is likely to 
be bo und. To estimate th e viri al mass they used M wu = 5Rcr 2 /G 
from Bertoldi & M cKed dl992l) . which we refined by deriving 
our own values for R = 0.3 pc (obtained from the interfe romet- 
ric map of dust emission from ICortes & Crutcherl l2006l as half 
the size of the continuum core or 10"), and cr = 1274 m s _1 
from our H I3 CO + observations (cr is not the FWHM value), ob- 
taining Mya - 563 Mq. Comparing this with the mass of the 
clump taken to be M submm = 3300 M© dCortes & Crutcherl2 006) 



6. Summary and conclusions 

We mapped the G30.79 FIR 10 molecular clump embedded in 
the W43 mini-starburst by mapping the molecular emission from 
12 CO(7 = 3 -» 2), I3 CO(7 = 3 -> 2), C I8 0(7 = 3 -» 2), 
CS(7 = 7 -> 6), HCO + (/ = 4 -> 3), H 13 CO + (7 = 4 -> 3), 
HCN(7 = 4 -» 3), and H 13 CN(7 = 4 -> 3). Emission from 12 CO 
is extended presenting a gradient from south-west to north-east, 
consistent with the dust morphology of the core. Even though 
we mapped 13 CO from a smaller region surrounding the cen- 
ter, its emission appears to follow I2 CO, which is also opti- 
cally thick. The C 18 emission was also mapped at the central 
30" x 30". Its line profile is almost Gaussian without any of 
the complex features seen in the previous CO isotopomers and 
is likely optically thin. We found a high velocity component at 
115 km s _1 , which was also present all over our sampled region. 
We interpreted this emission as a foreground cloudlet not asso- 
ciated with our clump. The emission from HCO + is optically 
thick, with an optical depth t = 11.3, but more compact than 
I2 CO. This is expected due to the higher densities required to ex- 
cite this molecule. Additionally, both CS and H 13 CO + are even 
more compact, with the emission coming mostly from the cen- 
tral 30" x 30" and consistent with the position of the dust max- 
ima obtained by others. By using our observations, we estimated 
the HCO + abundance X(HCO + )=2.4 x 10~ 10 , a depletion factor 
fn — 8 .4, which is consistent wi th estimations done toward other 
high-mass star forming regions dPurcell et al.ll2006l) . 

Outflows were looked for in the line-wings of the observed 
molecular lines toward the center of G30.79 FIR 10. All HCO + , 
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HCN, and CS profiles present excess emission in their line 
wings. However, the bipolar nature of this emission is inconclu- 
sive at the 3cr level. Infalling motions were also looked for by 
studying the profile asymmetry of our molecular observations. 
The blue asymmetry was estimated by calculating the normal- 
ized velocity difference (see Table [3}. The clearest evidence for 
infall is given by the HCN spectra. Toward the center and sur- 
roundings pointings, the emission is double peaked with the blue 
pea k stronger than the red peak (see Fig. 6 and fTOt. By using 
the Mvers et aD (1 19961) improved model by Di Francesco et al.l 
(120011) . we estimated an infall velocity of 0.5 kms 'with an in- 
fall rate of 5 x 10~ 3 M Q /yr. 

We refined a previous estimate for the magnetic field strength 
on the plane of sky in this region. Interf erometric observa- 
tions of polarized emission from dust from ICortes & Crutcherl 
(2006) estimated a magnetic field of B pos =1.7mG, by using the 
Chandrashekar-Fermi technique. The estimation needs a reliable 
tracer for the gas turbulent motions from the region traced by the 
dust emission. By using the line-width from our ion optically 
thin molecule H 13 CO + (7 = 4 — > 3), AV = 3 kms -1 , we were 
able to obtained an improved estimation of B pos =855 fiG for the 
magnetic field, which also refined the mass-to-magnetic flux ra- 
tio to A-l.9 or super-critical. Along with A, we calculated the 
contribution from turbulence to the virial mass, getting a value 
of Msubmm/Mvir = 5.9. These two results suggest that neither the 
magnetic field nor the turbulence have enough energy to coun- 
terbalance gravity. Therefore, the G30.79 FIR 10 clump must be 
bound and undergoing gravitational collapse. This also reinforce 
the core as an infall candidate. 
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